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Abstract 

One of the important approaches to detect quantum entanglement is using hnear en- 
tanglement witnesses (EWs). In this paper, by determining the envelope of the boundary 
hyper-planes defined by a family of linear EWs, a set of powerful nonlinear optimal EWs 
is manipulated. These EWs enable us to detect some three qubits bound MUB (mutu- 
ally unbiased bases) diagonal entangled states, i.e., the PPT (positive partial transpose) 
entangled states. Also, in some particular cases, the introduced nonlinear optimal EWs 
are powerful enough to separate the bound entangled regions from the separable ones. 
Finally, we present numerical examples to demonstrate the practical accessibility of this 
approach. 

Keywords rnonlinear optimal entanglement witnesses, mutually unbiased 
bases, MUB diagonal states 
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In the recent years it became clear that quantum entanglement [T] is one of the most important 
resources in the rapidly expanding field of quantum information processing, with remarkable 
applications such as quantum parallelism [2] , quantum cryptography [3] , quantum teleportation 
[H [5], quantum dense coding [6l [7] and reduction of communication complexity [8J. The 
above ideas are based on the fact that quantum entanglement, in particular, the occasionally 
occurrence of entangled states produce nonclassical phenomena. Therefore, specifying that a 
particular quantum state is entangled or separable is important because if the quantum state 
be separable then its statistic properties can be explained entirely by classical statistics. 

In this paper, we will deal with three qubit systems with 2'^-dimensional Hilbert space 
7^2 ® ® '^2, {'Hd denotes the Hilbert space with dimension d). A density matrix p on this 
Hilbert space, is called fully separable if it can be written as a convex combination of pure 
product states as follows 

p = j:p.|a«)(a«|«|af))(af|«|af)(af|, (1.1) 

i 

where la,-'''') are arbitrary but normalized vectors lie in 7^2, and Pi > satisfy YliPi ~ ^ 
(hereafter we will refer to fully separable states as separable ones). The first and most widely 
used related criterion for distinguishing entangled states from separable ones, is the Positive 
Partial Transpose (PPT) criterion, introduced by Peres [9]. Furthermore, the necessary and 
sufficient condition for separability in 0.2 ® Ti.2 and H2 ® TYs was shown by Horodecki in 
Ref. [To], which was based on a previous work by Woronowicz [11]. Partial transpose means 
transposition with respect to one of the subsystems. For a quantum state p^^ with matrix 
entries p^" = {ij\pAB\^n), the partial transposition with respect to the subsystem B, denoted 
by p^^, is defined by 

However, as it was shown in Ref. p^, in higher dimensions, there are PPT states that are 
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nonetheless entangled. These states are called PPT entangled states (PPTES) or bound 
entangled states because they possess the peculiar property that no entanglement can be 
distilled from them by local operations [T3]. Another approach to distinguish separable states 
from entangled ones involves the so called entanglement witness (EW) [14J. An EW for a 
given entangled state p is an observable W whose expectation value over all separable states 
is nonnegative, but strictly negative on p. There is a correspondence between EWs and linear 
positive (but not completely positive) maps via Jamiolkowski isomorphism [TB] . As an example 
the partial transposition is a positive map (PM). 

In this work, we consider those density matrices which are written as a linear combina- 
tion of maximally commuting observables taken from the set of tensor products ® Cg) 
C'', where A,B,C G {/2, ctx, o"?/, cTz} and i,j,k G {0,1} (cr^, ay and az are usual Pauh 
matrices). We will see later on that common eigenvectros of these observables form mu- 
tually unbiased bases {MUB)[16\ and so we will refer to a set of such observables as set 
of MUB observables, for instance by using the notation aicrjak = ai ® (Jj ® Ck, the set 
{III,az<7zI,<^zI<^z, I<^z<^z,<7x<7x<7x,<^x<^yO'y,o'yO'xO'y,o'yO'yO'x} IS a sct of MUB obscrvablcs. In 
fact, we consider tripartite MUB diagonal density matrices which are written in terms of 
MUB observables in a diagonal form. Then, we impose the PPT conditions (positivity of 
partial transposition with respect to all subsystems) to these density matrices and refer to the 
region of those density matrices which satisfy all of the obtained PPT conditions as "feasible 
region" (see Fig.l and Fig. 2 for example). In this way, we see that partial transposition plays 
an important role because in this type of density matrices, conditions obtained from positivity 
of partial transpositions are linear and feasible regions are completely contained in polygons; 
this allows us to investigate the separability or entanglement of the density matrices. In order 
to distinguish PPT entangled states {PPTES) from separable ones we construct some linear 
and nonlinear EWs. Namely, we consider density matrices that their common eigenvectors 
are maximally entangled states (GHZ-states) and construct an EW that detects such density 
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matrices. Finally, we consider three categories relevant to some special choices of the param- 
eters of density matrices, and by using the linear EWs wc distinguish the region of PPTES 
and separable states completely. In other words, for density matrices contained in one of these 
three categories, we show that if the introduced linear EWs can not detect their entanglement, 
then they are necessarily separable. 

We have also provided some numerical evidence suggesting that PPTES of three qubits 
can be detected by using the nonhnear EWs. 

The paper is organized as follows: In section 2, we introduce the MUB-{zzz)g diagonal 
density matrices and consider the corresponding PPT conditions and feasible region. Section 
3 is devoted to definition of an EW and construction of optimal hnear EWs. In section 4, 
we obtain an envelope of family of linear EWs and construct some nonlinear EWs. Section 
5 is devoted to classification and detection of PPTES for MUB diagonal density matrices in 
three categories. In section 6, we discuss some numerical analysis for evahiating the feasible 
region and the region of PPTES. The paper is ended with a brief conclusion together with two 
appendices. 

2 MUB diagonal density matrices 

In this section we introduce the so called MUB diagonal density matrices. The basic notions 
and definitions of MUB states relevant to our study are given in the Appendix /. 

2.1 MUB-{zzz)g diagonal density matrices 

In this subsection we introduce the MUB-{zzz)g diagonal density matrices which are consid- 
ered through the paper. A MUB-{zzz)g diagonal density matrix for three qubits is defined as 



8 



8 




(2.2) 



i=l 



1=1 
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where three qubit GHZ states {ipi) for i = 1, 2, 8 are given by 

1111)], 

1110)], 
1101)], 

|100)]. (2.3) 

Then, by using the MUB states in hne 6 of Table / given in Appendix A, the density matrix 
p can be rewritten as foUows 

p = + ri<7^aj + r2aJaz + r3lazCrz+r4^ax(T^ax+r5a^ayay + rGayaxCry + r7ayaya^], (2.4) 

o 

where 

ri = +Pl +P2+P3+P4-P5~P6-P7 - Ps, 

r2 = +Pi +P2-P3-P4 + P5+P6-P7 - Ps, 

r3 = +Pl +P2-P3-P4-P5-P6+P7+ PS, 

r4 = +pi -P2+P3-P4+P5-P6+P7 - Ps, 

r5 = -Pl +P2+P3-P4+P5-P6~P7+ Ps, 

Tq = -pi +P2+P3-P4-P5+P6+P7 - Ps, 

r7 = ~Pl +P2-P3+Pi+P5-P6+P7 - Ps- (2.5) 

It should be noticed that the MUB states of any line of Table / except for the states in the 
first three lines which are associated with separable states, can define a MUB diagonal density 
matrix, similarly. In the next subsection we impose the PPT conditions to the density matrix 
fl2.4p in order to obtain the corresponding feasible region (region of those MUB diagonal density 
matrices which satisfy all of the PPT conditions). 



1^1) 
iV'a) 

|V^5) 

iV^r) 



1 

V2 
1 

V2 
1 

71 



[|ooo) + |iii)], 1^2) 
[|ooi) + |iio)], 1^4) 

[|010) + |101)], 



[|ooo) 



-^[|011) + |100)], iV's) 



V2 

1 
1 



[|ooi) 



[|oio) 
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2.2 Feasible regions 



Here we are concerned with MUB diagonal density matrices of type (12 ■4p and by imposing 
the conditions obtained from positivity of partial transpositions with respect to each qubit, 
we obtain the so called feasible region. For these particular density matrices, the positivity of 
partial transpositions gives linear constraints on the parameters pi for z = 1, 2, 8. In order to 
obtain the feasible region for density matrix (12. 4p . first we group the conditions obtained from 
positivity of partial transpositions with respect to each subsystem (each qubit) in six partitions 

(P3,P4,P5,P6), {Pi,P2,P7,Ps), (Pi,P2,P5,P6), {P3, P4, Pi, Ps) , {pi, P2, P3, Pa) and {p5, Pg, Pi , Ps) as 
follows: 

The positivity of partial transposition with respect to the first qubit gives the following 
constraints: 

P3+P4+P5-P6>0 

P3+Pi-P5+P6>0 
P3-P4+P5+P6>0 
-P3 + P4 + P5 + P6 > 



(P3,P4,P5,P6) 



(2.6) 



{Pl,P2,P7,Ps) = < 



P1+P2+P7 -P8>0 
Pl+P2-P7+P8>0 

Pl-P2+P7+Ps>0 
-Pl +P2+P7+P8>0 

The positivity of partial transposition with respect to the second qubit gives: 



(2.7) 



(Pl,P2,P5,P6) = < 



Pl+P2+P5-P6>0 
Pl + P2 - P5 + P6 > 

Pl-P2+P5+P6>0 
-Pl + P2 + P5 + Pa > 



(2. 
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P3+P4+P7 -P8>0 

P3+Pi-P7+Ps>0 
P3-P4+P7+P8>0 
-P3+P4+P7+P8>0 

The positivity of partial transposition with respect to the third qubit gives: 



(2.9) 



{Pl,P2,P3,Pi) = < 



Pi + P2 + P3 - P4 > 
Pl+P2-P3+P4>0 
Pl-P2+P3+P4>0 
-Pl +P2+P3+P4>0 



(2.10) 



iP5,P6,P7,P8) 



(2.11) 



P5+Pfi+P7 -P8>0 
P5+Pfi-P7+P8>0 
P5-P6+P7+P8>0 
-Pb+Pfi+P7+P8>0 

The region of those density matrices of type (12. 4p which satisfy the above 24 constraints, is the 
feasible region. In order to specify the new perspective from this feasible region, we consider 
the parameters pi in four pairs (pi,P2), (^3,^4)5 (ps^Pe) and {p7,P8)- 

Now if we choose one of the pairs, say {pi,P2), then we can specify the projection of the 
feasible region to (^1,^2) plane with the following three inequalities ( the last inequalities of 
([23D, and (l210D . respectively PPT conditions) 

Pl<P2+P7+ P8 

Pi <P2+P5+ Pa 
Pi<P2+P3+ Pa 

By adding right hand side and left hand side of the above inequalities and using the equality 
Ylt=iPi ~ S^^ following inequality 



4pi - 2p2 < 1. 



(2.12) 
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Similarly, one can obtain the inequality 



Ap2 - 2pi < 1 



(2.13) 



by exchanging 1 and 2 in the above steps. For an illustration see Fig. 1. It should be noticed 
that if we chose any other pair from (^3,^4), (p5,P6) and (pjjPs) instead of the pair (^1,^2), 
we would obtain the similar inequalities as in (12.121) and (I2.13P for each pair. 
According to Fig. 1, since the vertex points (i, |), (^, 0), (0, j) and (0,0) satisfy the criterion 
(I2.12p . all of the points inside the shape will fulfill the PPT conditions (this is due to the fact 
that the feasible region is a convex region). 

We can find another new projection of the feasible region in the (^1,^3) plane, concerning 
the following inequalities 



I P3<Pa + P7+P8 

This region is illustrated in Fig. 2; therefore we have presented a projection of the spatial 
shape in a two-dimensional space. 

2.3 A special case of feasible region 

Here, we discuss a special case of feasible region which will be appeared in subsection 5.3 as 
a region of bound entangled MUB diagonal density matrices (PPT entangled states). To this 
aim, we consider the line pi +^3 = | of the feasible region (12.141) in the (^1,^3) plane (see Fig. 



First we take the feasible region for the (^3,^4) plane (see Fig. 3). If, we consider the 
following parametric line equation 



then (according to equations similar to ( 12. 12^ and ( 12.130 for the pair (^3,^4)), we obtain 




(2.14) 



2). 



1 

P3 = api + - 



(2.15) 



4p4 - 2p3 = 1 ^ P4 = y + I 



(2.16) 
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By substituting fl2.16p in f l2.15p and using the fact that < < 1/2 (see Eq. (12.141) ). one can 
obtain 

ap3 a + 1 a + 1 1 

Ps = => < Pa = < -, 

2 4 --^^ 4 -2a - 2' 

so we get — 1 < a < According to the boundary condition Pi+Ps = \ and (I2.15P we obtain: 

Pi = -api + ^. (2.17) 
Now by considering the PPT conditions 

~P3 +P4+P5+P6>0 , , 

(2.18) 

-P3 +Pi+P7+P8>0 

from equations (12.60 and (12. 9p . and adding the sides of the them, we obtain 

(1 - 2a)p4 - P2 > 0. (2.19) 

Also from the PPT conditions 

-Pl +P2+P5+P(i>0 
-Pl +P2+P7+P8>0 



(2.20) 



given in (12. 7p and (12.80 . we get 

- (1 -2a)p4+P2 > 0. (2.21) 

Then, from (12.190 and (I2.2ip we conclude that p2 = {I — 2a)pi. Then, by using (12.150 and 
(I2.17P one can easily conclude the following equations 

P3-P4 = {a- 1)P4 + ^ 

1 1 
Pi-P2 = -aPi + - - (1 - 2a)pi = {a- l)pi + -, 

which indicate that pz — Pi = pi — P2- On the other hand, from the inequalities (I2.19P and 
(I2.2ip . one can deduce that the left hand sides of the inequalities (12.180 and (I2.20p must be 
equal to zero. Therefore, for PPT density matrices with positive pi's, we obtain 

P5+P6=P7 + P8=P3-P4 = {<y-^)P4+^ >0 ^ Pa < 



4 - 4(1 - a) 
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furthermore we obtain > p^. Clearly, the PPT conditions f l2.10p and (12. lip are satisfied 
by using the relations Pi + p^ = P2 + Ps and P5 + Pe = Pv + Ps- The bound entanglement or 
separability of density matrices belonging to this special case of feasible region will be discussed 
in subsection 5.3 (as third category). For an illustration see Fig.3. 

2.4 'M.U'B-{zzz)g diagonal density matrices for which PPT condi- 
tions are necessary and sufficient for separabiUty 

In this section we consider the family of M\J'B-{zzz)g diagonal density matrices where the 
PPT criterions are necessary and sufficient for their separability. 

2.4.1 Case (1) 

In this case, we will put one of the pairs (^1,^2), (^3,^4), (psjPe) and {p7,ps) equal to (0,0), 
then we will see that the PPT conditions given in (I2.6I) - (I2.11I) are necessary and sufficient for 
separability of MUB-{zzz)g diagonal density matrices. For example, if we choose 

Pl=P2 = 0, 

then, by using the PPT conditions, we obtain 

P3 = P4, P5 = P6, Pi = P8- 

Then, the density matrices satisfying these conditions can be written as 

P = III 1 1 + na^aj + r2aja^ + rsla^cx^]. (2.22) 
o 

The density matrix p in (12.221) is separable, since by using (12. 5p we can rewrite p as 

P = ^{Psilll + (^zcrj - (TzIcTz - Icrz(^z) 
+p^{III + 0-2 Jcr^ - a^aj - I(TzO-z) + Prim + Io-zCTzI - ctzIctz - ctzCTzI)} = 

P3(IV'3)(^3| + \H{H)+P5{\^5){H + \^e){M+P7i\i^7){H + \i^s){H), 

which is clearly a separable state, since it is a convex combination of projection operators. 
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2.4.2 Case (2) 

In this case, we choose p^'s in each pair except for one of them to be equal, then we show that 
the PPT conditions (I2.6p - fl2.1ip are necessary and sufficient for separabihty of MUB-{zzz)g 
diagonal density matrices. For example, we consider 

Pi 7^ P2, P3 = P4, Pb = P6, Pi = Ps- 
Then, we can write the density matrix ( 12. 2p as follows 

pd\^B){H + \A){A\)+P7{\^7){H + \^s){H)- (2.23) 
We assume that Ps < P5 < P7 and p2 < pi (the other cases give the same results as those which 
is obtained by this assumption in the following). 

By substituting = in PPT conditions fl2.10p . we obtain pi < P2 + 2^3, so that we can 
write pi = P2 + 2p3 — 2ei (pi > P2 ^ < ei < pa). Also, from the assumptions ps < ps and 
Ps < P7, one can write ps = P3 + 65 and pi = ps + e-j, respectively. By substituting these values 
of Pi's in the density matrix (12.231) and using the resolution of identity Ylt=i = ^^^^ 

one can write 

p = eJII + (p2 - ei)(|^i)(^i| + \iJ2){H) + iP3 - ei)(/// + - |V^2)(^2|) 

+ e5(|^5)(^5| + |^6)(^6|) +e7(|^7)(^7| + |^8)(^8|)- (2.24) 

Then, from the fact that (/// ± \ipi){4'i \ — \4'2){4'2\) are separable states, one can see that for 
ei < p2 in f l2.24p . the density matrix p is separable (since it is written as a convex combination 
of product states). For ei > p2, we can write p as follows 

p = P2///+ (ei -p2)(///- - \^p2){H) + iP3 - ei)(///+ - |^2)(V'2|) 

+e5(|^5)(^5| + |^6)(^6|) + e7(|^7)(^7| + l^^s) (H) , 

which is again a separable state. So in the second case, the PPT conditions are necessary and 
sufficient for separability, too. 
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3 Entanglement witnesses 

An entanglement witness acting on the Hilbert space Ti. = 7^2 ® ® 7^2 is a Hermitian 
operator W= , that satisfies Tr{Wps) > for any separable state ps in BiTi) (Hilbert 
space of bounded operators), and has at least one negative eigenvalue. If a density matrix p 
satisfies Tr(Wp) < 0, then p is an entangled state and we say that detects entanglement 
of the density matrix p. Note that in the aforementioned definition of EWs, we are not worry 
about the kind of entanglement of the quantum state and we are rather looking for EWs which 
possess nonnegative expectation values over all separable states despite of the fact that they 
possess some negative eigenvalues. The existence oi an EW for any entangled state is a direct 
consequence of Hahn-Banach theorem [20] and the fact that the subspace of separable density 
operators is convex and closed [21]. Geometrically, EW^ can be viewed as hyper planes which 
separate some entangled states from the set of separable states and hyper plane indicated as 
a line corresponds to the state with Tr[iyp] = 0. 

Based on the notion of partial transpose, the EWs are classified into two classes: decompos- 
able EW^ {d-EW) and non-decomposable EW^ {nd-EW). An EW W is called decomposable 
if there exist positive operators P, Qk so that 

W = P + Ql^ +Ql^ +Ql% (3.25) 

where Tk, K = A,B,C denotes the partial transposition with respect to subsystems A, B 
and C, respectively. W is called non-decomposable if it can not be written in this form [22] . 
Clearly a d-EW can not detect bound entangled states (entangled states with positive partial 
transpose (PPT) with respect to all subsystems) whereas there are some bound entangled 
states which can be detected by a nd-EW. 

Usually one is interested in finding optimal EWs W which detect entangled states in an 
optimal way. An EW W is said to be optimal, if for all positive operators P and e > 0, the 



nonlinear Entanglement Witness using MUB 



14 



new Hermitian operator 

W' = {l + e)W-eP (3.26) 

is not anymore an EW [23]. Suppose that there is a positive operator P and e > such that 
W = (1 + e)W — eP is yet an EW {Tr{W'ps) > for all separable states ps)- This means 
that if Tr{Wps) = 0, then Tr{Pps) = 0, for all separable states ps which indicates that, the 
operator P is necessarily orthogonal to the kernel of W denoted by Ker{W). By using the 
fact that every separable state is convex combination of pure product states, one can take ps 
as a pure product state Also, one can assume that the positive operator P is a pure 

projection operator, since an arbitrary positive operator can be written as convex combination 
of pure projection operators with positive coefficients. 

3.1 EWs detecting bound MUB diagonal density matrices 

By employing tensor products of pauli operators relevant to M\JB-{zzz)g state of Table / of 
the Appendix A, we introduce the following linear three qubit EW pi] 

W = AqIII + Ai la^a^ + y42(cr^cr^a^. + a^-o-yO-j^) + 743(0-^0-^.0-^^ + o-,yayO-^.), (3.27) 

where 740,^2,^3 > and Ai can be negative or positive. Now evaluating the trace of EW 
f l3.27p over a pure product state, 

Ps = \a){a\0\P){P\0\-f){j\, 

we get 

Tr[Wps] = Ao + A163C3 + A2{aibiCi + aib2C2) + ^3(026102 + a262Ci), 

where aibjCk := Tr[ai(jjakPs\ for i,i,k = 0,1,2,3 with oq = &o = cq = 1. We parameterize 
points on the unit sphere S'^ using traditional spherical coordinates, so that 9 and (p stand for 
the angles of colatitude and longitude, respectively {9 G [0,7r],y9 G [0,27r]). Thus, the points 
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a = (01,02,03), b = (&i,&2,&3) and c = (01,02,03), can be uniquely represented as the unit 
vectors with the following coordinates 

oi = sin6'i cosyji, 02 = sin6'i siny^i, 03 = cos6'i 

bi = sin 02 cos ip2, ^2 = sin 6*2 sin ip2, &3 = cos 62 
Ci = sin 63 cos ips, C2 = sin 6*3 sin ^93, 03 = cos 6^3, 

so that, we obtain 

Tr[iyp5] = Ao+Ai cos 6^3 cos6'2+sin6'i sin 6*2 sin 6*3(^2 cos(/Ji cos {(p2 — V'3)+^3 sinyji sin {ip2 + V^3))- 
By appropriate choice of the angles, one can minimize the above expression. In fact, in the 



Appendix B, it has been proved that by taking Aq = + and Ai = —y^A^ + Al, the 
minimum value of Tr[l^ps] is attained to zero, i.e., we obtain 

min(Tr[W ps]) = 0. 

Consequently, EW fl3.27p takes the following form 



VT^'^^ = y Al + Al{III - Ia;,az + cos'ilj{a^a^a^ + a^ayay) + smip{aya^ay + ayayax)), (3.28) 
where 

, A2 . , A3 

cos?/; = — , smip 



In the following, we discuss the optimality of the obtained linear EW W^. 
3.2 Optimality of the linear EW VF^'^) 

According to the arguments about optimal EWs given in section 3, in order to prove the 
optimahty of the EW W^'f'^ given in (13.281) . it suffices to show that there exists no positive 
operator P such that W := (1 + e)W'^'^^ — eP be an EW, namely it must be proved that for 
any pure product state |z/) so that Tr{W'^''^^\v){i'\) = 0, there exists no positive operator P 
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with the constraint Tr(P|z/) (z/|) = 0. By considering a general three qubit pure product state 



as 



one can evaluate 



V2 



cos [^j 



cos(f) 
eiv2 sin { §2 



( 



2) J 



cos(f 



oi^i sin 



(3.29) 



Tr\W'^'^^ Ps] = 1— cos 02 cos6'3— sin 6i sin 6*2 sin 6'3(cos^/' cos ipi cos {ip2 — v53)+sin ip sin ipi sin {ip2 + v^a)), 
where ps = Now, it is easily seen that by choosing the angles 9 and 9? as follows 

cos (<^2 - (^S) = 1 TT ^ ^ TT ^ ^ 

sin {(p2 + (ps) = 1 



(2): 



cos (v92 - v^a) = 1 
sin {ip2 + (ps) = -1 



TT TT 



COs((^2-<^3) = -1 TT 37r / ^ ^ 

(3) : <( ^ ¥?2 = ^, ¥^3 = - V'l = ^ - ^> ^1 = 2 ' ^2 = ^3, 

sin (v32 + V'3) = -1 



(4): 



C0s{ip2 - ips) = -1 



StT TT TT 

V52 = -^^3 = --,(pi = n - ip,ei = -,02 = 03, (3.30) 



sin {(p2 + (p3) = 1 

we obtain TrlW^"^^ ps] = 0. Now, in order to prove that W^'^^ is an optimal EW, we proceed 
as follows: Let P be a pure projection operator that one can subtract from W^'^\ so that 
(1 + e)W^'f'^ - eP is an EW for some e > 0. From Eq. f l3.28p . one can easily see that, any 
pure state of the form = \a)\z^z^) + \(3)\z_z„), (where, |a) and \P) are arbitrary states) 
belongs to the KeriW^"^^), i.e., we have Tr[iy('^)|$)(^|] = 0. Then, due to the fact that, the 
pure projection operator P must be orthogonal to Ker{W^'^^) (and so orthogonal to |^)(^|), 
we have P = |$)($| with 

= \a)\z^z,) + \(3)\z_z+) (3.31) 
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Now, the pure projection operator defined as above, must be orthogonal to pure product states 
li/*^*-*), i = 1,2,3,4 obtained by substituting the angles given by fl3.30p in fl3.29p . since these 
states belong to Ker{W^^^). But, this is possible only if ai, 02, Pi and P2 satisfy the following 
equations: 



(z/i|<l>) = 


(«! + 0:26" 


+ (A + 1326-''^) 


= 0, 


(1/2!$) = 


(ai + 0:26^ 




= 0, 


Z/3|$) = 


— (tti — 




1 = 




— (ai — 


+'^) + (A - /32e+^^; 


) = 



Above equations imply that for 

^ ^ 0,7r 

we have 

ai = a2 = Pi = f32 = 0. 

Therefore, there is no positive operator P to subtract from W^'^\ 
In general, linear optimal EWs can be written as 

Wi%j,kui,m) = III ± + cos^(0, ± Ok) + sin^(Oi ± O^), (3.32) 

where i = 1, 2, 3 while the indices j ^ k ^ I ^ m take values between 4, 5, 6, 7. 
The observables Oj for i = 1, 2, 3 and 0^ for j = 4, 5, 6, 7 are defined as 

Oi = la^a^, O2 = crja^, O3 = a^aj, 



4 Non-linear optimal EWs 

Actually with a given entangled density matrix, one can associate a non-linear EW, simply 
by defining a non-linear functional, so that it is nonnegative valued over all separable density 
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matrices, but it is negative valued over the density matrix. In other words, we optimize 
Tr[W^^^^^j^j^^^^i^^^p] where, are the hnear optimal EWs given by fl3.32p and p is 

the MUB-{zzz)g diagonal density matrix given by (12 ■4p . Then, one can easily get 

Tr[W^^^^^.^f^^^^i^^^p] = (1 ± Ti) + {vj + Tk) cosip + {ri + r.^) sinifj, 

which indicates that, by appropriate choice of the parameter ^ as a functional of p, one can 
obtain a non-linear function of the parameters of p which is nonnegative over all separable 
states. To this aim, we define 

n + 
cos 9 



then Tr[W^^^,. ^^ ^p] can be written as 



^/{rj + + {ri + r„)^ 

±^Mlml,m)p\ can be writtei 

Tr[l^(t,La,fc),(^n^)P] = 1 ± + sJ{r,+rkf + {ri + rmYcos{i^ - 9). 
Now, by choosing {ip — 6) = vr, we obtain 



Tr[W±i^±y^kUi,m)P] = l±ri- {r^ + r^f + {n + r^f . (4.33) 

The above expression is the required non-linear function in terms of the parameters of p and it 
is definitely nonnegative valued function of separable states, hence it is the non-linear optimal 
EW associated with p (since it is obtained from optimal linear EWs). 

4.1 Non-linear EWs as an envelop of family of linear EWs 

As the parameter ijj of linear EWs k) {i m) varies, the envelope of hyper planes defined 

by 

namely their intersections, define the boundary of PPT bound entangled states that can be 
detected by the linear EWs. Obviously, the envelope of these curves can be obtained simply 
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by eliminating the parameter ip from the Eg. (14.341) . To this aim, we need to determine cosip 
and sin ip by solving above equation together with the equation that can be obtained by taking 
its derivative with respect to ip equal to zero, i.e., we consider 

Tr[W^^^^^^^^^^^Qj^p] = (1 ± Ti) + (r4 + rs) cosip + (rg + Vj) simp = 0, 

ip{Tr[W^^^^^^^^^^^^^^^p]) = - sin^in + rg) + cos^(r6 + ry) = 0. 

By solving the above equations one can obtain 

, {vj + rk){l ±ri) 

[rj + rky + [ri + VmY 

. , +r„)(l ±ri) 

smip = — —. 

[rj + rkY + (n + rmY 

Now, using the identity cos^ ip + sin^ ■?/' = 1 we obtain the required envelope of curves defined 
by the following equations 

{l±rif = {r^+rkf + {ri + r^f. 



5 Bound entangled MUB diagonal density matrices 

In this section, we consider three main categories of bound entangled states according to 
equations (IB-ip and (IB-iip given in the Appendix B. In these categories, the relations \{rj ± 
rfc)| < (1 ± Tj), are always satisfied since if we consider for example the inequality 

\{r^-r^)\ > (1 + ri), 

then we conclude the inequality P2 + P4 < which is clearly impossible. 

5.1 First category 

The first interesting family of three qubit bound entangled states is introduced for the choices 
of the parameters so that: 

1 ± ri = rs ± rg, 1 ± ri = r^^ r^,! ± r2 = r^ ± ry. 
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1 ± r2 = r4 =F re, 1 ± rs = re ± ry, 1 ± rs = r4 =F r^, 



(5.35) 



for example, if we consider 1 + ri = — ry, then we obtain P2 = Pi = 0. Then, the PPT 
conditions (I2.6l) - fl2.11l) lead to pi = p^, and triangle inequalities for the cases (pi,P5,P6) and 
{PitPttPs) are estabhshed. This state can be detected by non-linear EW W^+i,-(4,7),(5,e)5 since 
by using the result (14. 331) . we have 



which indicates that for 7^ re, rr[VF+i _(4 7) (5 e)p] < 0. 

On the other hand, by imposing the condition (r5 — tq) = P5 — Pe — Pr + ps = 0, the state 
p will be separable, since we have 



So, by using the relations fl2.5p . one can see that if r2 > (r2 = 2(pe — Ps) since we have 
Pi = Ps and P2 = P4 = 0), then we get pe > Ps and therefore we have 

p = (1 - r2)/// + (ri + r2)I(T^az + r2[(/// - a^aj){in + (x^/a^)] + r^a^a^cfJx + r^a^ayay 

with 



r4 — ry = 4pi 
ri +r2 = 2(pi -pr-ps) 
ta-tt - (ri + r2) = Api - 2pi + 2p^ + 2p8 = 2(pi + + Ps) 
1 - r2 = 1 - 2p6 + 2p8 = 2(pi + P7 + Ps) 
1 - r2 - (ri + r2) + r4 - ry = 4(pi + pj + ps) < 1, 



which is separable state. 

If r2 < then we will have ps > Pe and 

p = (1 + r2)/// + (ri - r2)I(Tz(T^ - r2[{III - (J^aJ){in + oja^)] + r^a^a^a^ + r^a^aya. 



Tr[W^+i,_(4,7),(5,6)P] = (1 + ri) - v/(l + ri)2 + (r5-r6)2. 



(5.36) 



P5 - P6 - P7 + Ps = ^ P5 + Ps = P6 + P7- 
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with 

ri - r2 = 2{pi -p^-Pe) <0 

r^-r-r - {n - = 2(pi + ps + Pe) > 

l+r2 = 2(pi+p5 + P6) >0 

(1 + r2) + r4 - rr - (ri - r2) = 4(pi + ps + Po) < 1- 

For the above cases we had P2 = P4 = 0, so, this category consists two vanishing non-paired 
{p2 and p4 belong to different pairs) parameters. The other cases in fl5.35p can be discussed 
similarly. 

A special case 

As a special case, if we consider p in (12. 4p with the following parameters 

P4 = P8 = P6 = 0, 
P3=P5=P7 = P, 

we get Ti = r2 = and = = r^. Then, p can be written as 

p = + ri{Ia-^az + a^/cr^ + a^aj) + r^an^a^ax + r5{axayay + aya^ay + CTj^o-yO-x.)]- (5-37) 

Concerning the following PPT and normalization conditions 

Pi+P2-Pz>^ 

Pl-P2+P3>0 
-Pl +P2+P3>0 

8 

y^Pi = 1 ^ Pl + P2 + 3p3 = 1 

we construct the convex hull of the following boundary planes 
Pl + P2 - P3 = 
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Pi + P2 + 3p3 = 1 ^ 4p3 = 1 p3 = I 

Pl-P2+P3 = 

Pi + P2 + 3p3 = 1 =^ 2pi + 4p3 = 1 
- Pi + P2 + P3 = 

Pi + P2 + 3p3 = 1 ^ 2p2 + 4p3 = 1, 

which define a triangular bound entangled region in (pi,P2,P3) space (as it is shown in 
Fig. 4). With the boundary defined by the following fines where the boundaries of triangular 
region corresponding to the lines passing through the points 0, |), (0, i, i) and (i,i,0), 
where the states corresponding to the sides defined by the thick fine passing through the 
points {\,0,\) and (|, |, 0) are separable. 

5.2 Second category 

Another interesting family for bound entangled states is obtained by considering the following 
cases: 



l±ri 


= r4 + r5, 1 ± ri 


= ^4 


- r5, 1 ±ri 


= re + ry, 1 ± ri 


= r-Q 


- 


l±ri 


= r4 + r6,l±ri 


= r4 


- re, 1 ± ri 


= rg + r7, 1 ± ri 




- r7 


l±r2 


= r4 + r5,l±r2 


= r4 


- rs, 1 ± r2 


= re + r7, 1 ± r2 


= re 


- rj 


l±r2 


= r4 + ry, 1 ± r2 


= ^4 


- r7, 1 ± r2 


= rg + re, 1 ± r2 


= r5 


- re 


l±r3 


= r4 + re, 1 ± rs 


= r4 


- re, 1 ± r3 


= rg + r7, 1 ± r3 


= '^s 


- r7 


l±r3 


= r4 + r7, 1 ± r3 


= r4 


- r7, 1 ± r3 


= rs + re, 1 ± r3 


= r5 


- re 



If we choose one case such as 1 + ri = r4 + re, we obtain p? = Pi + P2 + 2p4 + pg. Then, by 
using PPT conditions we obtain 

P1+P2-P7+P8 >0=^P4 = 0, (5.39) 

P3-P7+P8 >0=^P3=Pl+P2, (5.40) 
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and satisfy the triangle inequality for (p3,P5,P6) case. So according to f l5.39p and (15.401) we get 

P7=P3+P8- (5.41) 
Applying the EW M/+i,+(4,6),(5,7) to the state ( 15.37p . we obtain 

Tr[Wp\ = (1 + ri) - V(l + ri)2 + (r5 + r7)2 < 

the condition {r^ + rj) = —pi + P2 + Pb — Pe = corresponds to separable state. 
By using the normalization condition 

8 

J^Pi = 1 ^ 3p3 + P5 + pe + 2p8 = 1 

i=l 

and 

-Pi + P2 + P5 - Pe = ^ pi + pe = P2 + P5, 

we obtain for r2 > 

p = (1 - r2)/// + (ri + r2)I(J^a^ + r2[(/// + aja^){in - a^aj)] + r4(Ta:Cra;Cr3; + r^aya^ay 

r4 + re = 4(pi + ps) = 4p3 

^1+^2 = -4j98 

r4 + re - (ri + r2) = Apj 
1 - r2 = 4p7 
^ (1 - r2) + (r4 + re) - (n + rs) = 8^7 

while for r2 < 0, we get 

p = (1 + rs)/// + (ri - r2)Icr^a^ - r2[{III + a^aJ){III - aja^)] + r^a^a^a^ + r^aya^ay 

ri - r2 = 2(p3 - P5 - Pe) < 
r4 + re - (n - r2) = 2(p3 + + Pe) > 
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l+r2 = 2(l-2p7) =2(P3+P5+P6) >0 

^ (1 + r2) + (r4 + re) - {n - = 4(p3 + P5 + pe) < 1. 

So, in this case according to (IB-ip . if we choose parameters as Pi = pj + pk, {pi and pj are 
in the same pairs) and pk belong to another pairs, then we obtain second category of bound 
entangled states. 

The other family (15.21) can be considered similarly. 

5.3 Third category 

The last category is given by the following cases: 

1 ± ri = r4 ± ry, 1 ± ri = rs =p rg, 1 ± r2 = r4 ± rg, 

1 ± r2 = rs =F ry, 1 ± rs = r4 ± rs, 1 ± rs = re =F r^. (5.42) 

If l-n = r4-r7 thenpi+p3 = P2+P4+P5+P6+P7+P8 = I and rs-re = 2(p5-p6-P7+P8)- 
The PPT conditions for this case have been previously considered (section 2.4.2). Therefore, if 
r^ — TQj^O this state will be bound entangled and can be detected by W^-i,-(4j),(5,6); otherwise 
it is separable since we have 

r5 - re = ^ Pe + Pt = P5 + ps 
and so from the PPT conditions we get 

P5 =P7,P6 =P8 

after calculation of rj's, we obtain 

p = 2p2{III + cx,aJ){III + I a, a,) + 2p4(/// + cx,aJ){III - a J a,) 
+ (1 - 2{p2 + Pa))! II + 2p^a^(y^(j^ - Ap^a^ayay. 
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We know that the first and second cases in the density matrix are separable and for other cases 

|(1-2(P2+P4))| + |2P5H-|-4P6|<1. 

So, if we choose two parameters and add them as Pi + Vj = \ ;( Vi^Vj are in different pairs) 
then, the category consists the bound entangled states. 
The other cases fl5.2p can be discussed similarly. 

6 Numerical analysis of three-qubit bound MUB diag- 
onal entangled states 

This section is devoted to some numerical studies of three-qubit bound MUB diagonal entan- 
gled states as follows: The feasible regions in (^1,^2) and (^1,^3) planes, defined by equations 
(]2.12p and ( I2.14p . are supported numerically. Using the nonlinear EWs, about 2.7% of bound 
M\J'B-{zzz)g diagonal density matrices are detected numerically. The numerical results are 
plotted in (^1,^3,^5) , (P2,P4,P8) , (Pe,^?) and (^1,^3) , (^2,^4) , (P5,P6) , {.Pi^Ps) phase spaces 
and the bound density matrix (pi = 0.043425,^2 = 0.15308, pa = 0.016132,^4 = 0.19387, = 
0.059793,^6 = 0.24806, p7 = 0.18207, pg = 0.10357) is shown in Fig.5 and Fig.6, as a prototype 
of a bound MUB diagonal density matrix. 

7 Conclusion 

The feasible region of PPT MUB-{zzz)g diagonal density matrices is determined, where it is 
a convex polyatope due to linearity of PPT conditions. In order to detect three-qubit bound 
MUB diagonal entangled states, some nonlinear optimal EWs are manipulated, such that they 
form the envelope of the boundary hyper-planes defined by a family of optimal linear EWs. 
By using these nonlinear EWs, the region of bound entangled states and separable ones are 
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determined analytically in some particular cases three categories, where the numerical analysis 
support them. The results thus obtained in this paper indicate that, the proposed methods 
in this work, can be used in studying entanglement of more general systems with linear PPT 
conditions such as multiqubit MUB diagonal systems. . 



Appendix A 



I. Mutually unbiased basis 

Let F be a (/-dimensional Hilbert space with two orthonormal basis 

5i = {|ei), |e2), |ed)} and 

52 = {|/l),|/2),...,|/d)}, 

where \ei) and \fi) for i — 1, 2, d belong to C*^ (the standard Hilbert space of dimension d 
endowed with usual inner product denoted by ( | )). 

The basis Bi and B2 are called mutually unbiased if and only if 

= (A-i) 

As an example, for a two-level system there is such a set of bases that can be represented in 
terms of eigenvectors of the usual Pauli matrices cr-r, Uj,, as follows 

5. = {^(|0) + |1)),^(|0)-|1))}, 

^. = {;^(|0)+^|1)),;^(|0)-^|1))} and 
5. = {|0>,|1>}. 

When c? is a prime or power of a prime, the maximum number of such MUB's is equal to d+ 1, 
otherwise there is no clear number of sets. 
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According to Refs. [HI [19] and Table /, in the case of three qubits, we have nine sets of 
mutually unbiased bases and corresponding maximally commuting sets of observables, where 
we will refer to them as generalized Pauli matrices; each of these sets consist of seven commut- 
ing observables. In the Table I, the first three rows contains product common eigenvectors, 
{xyz).„, (yzx)-,, and {zxyj^, (subscript n means product state). For example, eight states for basis 
{xyz)-„ could be written ) where = 1 and rij = correspond to spin down and 

spin up along the ith axis for i = x,y, z, respectively. These product states are separable so 
we will not use them for construction of EWs. Other bases consist of six maximally entangled 
states {xxx)Gi, {yyy)G, {zzz)g, {xzy)G, {yxz)G and {zyx)G, here subscript G denotes a family of 
Greenberger-Horne-Zeilinger {GHZ) states. For example, eight states for basis {zzz)g and 
{xxx)Gi can be written as 



where labels bar show that if = or 1, then = 1 or 0, respectively. In section 2, we have 
considered only the state {zzz)g', since the other cases can be obtained from this state by local 
unitary operations. 

II. MUB sets for three qubit systems 

We know MUB can be constructed using a number of methods that depend on the dimension- 
ality of the space. These methods using for different case such as dimension space is prime, a 
product of primes, or a power of a prime, and if it is odd or even. We confine our study to the 
case of three qubits, that is, to an eight-dimensional Hilbert space, in this space there exist 
four MUB formation, where denotes sets of MUBs where the basis vectors are either separable 
, biseparable or entangled states. The four MUB formations are (2,3,4) (here 2 means two 
separable states, 3 means there is three biseparable state, 4 means four maximally entangled 
states), (0,9,0), (3,0,6) and (1,6,2). There are other MUB sets for three qubits but one can 



{zzz)g = \n^n^n„ ±) = {\n,n,n,) ± \n,n,n,)), n, = 0,l 



(A-ii) 



{xxx)Gi = \n^n^n^, ±) = {\n^n^n^) ± i\n^n^n^)), = 0, 1 



(A-iii) 
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get them by local transformations from the previous ones. In this paper we work with table 
(3, 0, 6), line (6) in order to introduce our EWs. This table is given by 



1 


(xyz)^ 




lOyl 


IIcTz 


a^ayffz 


a^CTyl 


(^xIcTz 


IcXyaz 


2 


(yzx)^ 


ayll 


lOzI 


n^Tx 


(Tyazax 


(TyCFzI 


aylcTx 




3 


(zxy)^ 




lOxl 


Hay 


(yz(yx<yy 


(^zC^xI 


a Jay 


la^CTy 


4 


{xxx)Gi 




OzOyOz 


a^a^^Gy 


ayUyGy 


(^x<^xl 


(^xlcyx 


I<yx<^x 


5 


{yyy)G 




(^xC^zC^x 


(7x(7x(7z 


(7z(7z<^z 


(Tydyl 


(Tyldy 


Idydy 


6 


{zzz)g 


OxOyOy 


(yy(yx(yy 


(yy(yy(yx 


(^x<^x<^x 


O-zCTzI 


(Tzlf^z 


la^cTz 


7 


ixzy)G 


(yzC^xOz 


(Ty(Tx(Tx 


(yy(Ty(Tz 


(^z(^y(^x 




a^IcTy 


la^CTy 


8 


{yxz)G 


OxOyCTx 


(TzCTyCTy 


(^z<^z<^x 


(^xCTzCTy 


aya^I 


CTylaz 


Icrx<^z 


9 


{.zyx)G 


OyO-^Oy 


(^x<^z<^z 


(^x<^x<yy 


GyG^Gz 


a^Gyl 




laya^ 



Table 1: Nine sets of operators defining a (3,0,6) MUB . 
In this table three states 



(A-iv) 



{xxx)Gi, iyyy)G, {zzz)g 

can be reversibly converted into each other by local unitary operations (permutation), (i.e. 
Cy ^ ^ (Tx), and e.g., if we construct EW using {zzz)g then this EW can be converted by 
another EW for {xxx)g state by applying the local unitary transformation . 
Another states 

{xzy)G, iyxz)G, izyx)G (A-v) 
can be transformed into each other by local unitary operations, e.g., for state {xzy)G we have 



{xzy)G = \nlnlnl,±) = (|ra>^^J) ± \nlnlnl)), 



(A-vi) 



we can convert the states flA-vP into other three maximally entangled states flA-ivP by the 
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following local operations 



U = ^ 



^ 1 ^ 



v 



J 



\ 



ex 
e-x 



U = ^ 



1 I 



where, we have applied the permutation x ^ z to the first qubit, y ^ x to the middle qubit 
and y ^ z to the rightmost qubit. 

The above discussion was used for table (3,0,6), but if we choose table (2,3,4), we will 
obtain maximally entangled states {{xzy)Q, {yyz)G, {yxy)G, {zyx)Gi) so that, according to the 
following local operators and permutations we can convert the states (\A-v\i into other three 
maximally entangled states. 

{yyz)G {yxy)G with local operators [/(l)^^^, t/(2)y^^, f/(3)^^y .(the notation U{i) 
means that U acting on the (i)-th qubit) 

{xzy)G — ^ {zyx)Gi with permutation {ax — > — > (Ty) only for first and second qubits and 
local operators U {S)y^x 

{xzy)G — * {yyz)G with local operators U{l)x^y, U {2)z^y and with permutation [ay — — ^ 
a^) for the rightmost qubit. 



Appendix B 



The following cases are used for construction of EWs: 



r4 + rs = 2(j93 -PA + Pb- P&),rQ + r-j = 2{-pi + P2 + Pj 



(B-i) 



r4 - rs = 2{pi -P2+P7- P8),rG - rj = 2{p^ - Pi - P5 + Pg) 
r4 + re = 2(p3 -p^ + pj- Pa), + rj = 2{-pi + P2 + Ps - Pe) 
r4 - rg = 2{pi -P2+P5- P6),r5 - rj = 2{p^ - Pi - Pr + Ps) 



r4 + r7 = 2{p5 -P6+P7 - P8),r5 + = 2{-pi +P2+P?,- Pi) 
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n-r-j^ 2(pi - P2 + P3 - P4), - re 
We choose one the following cases for our EW: 

1 + ri = 2(pi +P2 +P3 +P4), 1-ri 

1 + r2 = 2(pi + P2 + P5 + Pe), l-r2 
1 + rs = 2(pi H-p2 + P7 +P8), 1 - rs 
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2(P5 -P6 -P7+P8)■ 
2(P5+P6+P7+P8) (B-ii) 

2(P3+P4+P7+P8) 
2(P3 +P4 +P5 +P6)- 
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Proof for EWs detecting bound 'M.\J'B-{zzz)g diagonal density matrices 

For more detail, we have 

Tr[iyps] = Ao+Ai cos 6*3 cos 6*2+8111 9i sin 02 sin 93{A2 cosipi cos {(p2 — ^Jsj+As sin ipi sin (932 + ^s))- 
Taking = 9^3 = t , 6'i = f 



and according to relation {—^/a^ + < asin^ + bcos9 < y/a^ + IP') , we have 

Tr[Wps] ^ Aq + Ai cos 6*3 cos $2 ± sin 02 sin Os^A^ + Al, 
according to cos ^3, we obtain 

Tr[Wps] = A ± ^A?cose2' + sin^2'(^i + ^i), 
and from the condition (A^ + A| > A^, ^2 = f ), one can get 



Tr[Wps] = Ao ± ^JAl + Al 

so , if 



Ao = ^JAI + Al 
Tr[Wps] > 0, 
W ^^Al + AlIII + Ai Ia,a,-A2{ 
then Al is an arbitrary number, therefore by taking Ai = — A3^^ A3, we get 



+ ^3 V + M 
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Figure Captions 

Figure-1: Shows the feasible region in (^1,^2) plane defined by Eq. fl2.12p . 

Figure-2: Shows the feasible region in ( ^1,^3) plane defined by Eq. fl2.14p . 

Figure-3: The two colored triangles show the feasible region obtained by the condition 
Pi + Ps = where the light gray-colored triangle determines the region in (^3,^4) plane and 
(according to the PPT conditions) the dark gray-colored triangle indicates the region in ( 
Pi,P2) plane. The dotted line represents parametric equation (12.151) which has been drawn for 
two regions, i.e, from PPT conditions, dotted line in ( ^3,^4) plane corresponds to the dotted 
line in ( pi,P2) plane. 

Figure-4: The gray-colored triangle shows the region of bound entangled states in {pi,p2, Pa) 
phase space and thick line shows the region of separable states. 

Figure-5:(a): Perspective of the region of bound entangled states in (^1,^3) plane obtained 
by numerical analysis, (b): Perspective of the region of bound entangled states in (^2,^4) 
plane obtained by numerical results, (c): Perspective of the region of bound entangled states 
in {p5,Pq) plane obtained by numerical analysis, (d): Perspective of the region of bound 
entangled states in {p7,p%) plane obtained by numerical analysis. 

Figure-6:.(a): Perspective of the region of bound entangled states in (^1,^3,^5) phase 
space obtained by numerical analysis, (b): Perspective of the region of bound entangled states 
in {p2,P4:,P8) phase space obtained by numerical analysis, (c): Perspective of the region of 
bound entangled states in (pgjP?) plane obtained by numerical analysis. 



